ITP-SB-95-16 



Bottom quark production cross section at 
fixed-target pp experiments 



N. Kidonakis and J. Smith 

Institute for Theoretical Physics, State University of New York at Stony 
Brook, Stony Brook, NY 11794-3840, USA 

June 5, 1995 
Abstract 

The cross section for bottom quark production at fixed-target ener- 
gies is calculated for a wide range of beam momenta. A detailed anal- 
ysis is given for the HERA-B experiment. We consider both the order 
Oig cross section and the resummation of soft gluon corrections in all 
orders of QCD perturbation theory. The inclusive transverse momen- 
tum and rapidity distributions, including resummation, for bottom 
quark production at HERA-B are also presented. 
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1 Introduction 



The calculation of production cross sections for heavy particles in QCD is 
made by invoking the factorization theorem [|IJ and expanding the contribu- 
tions to the amplitude in powers of the coupling constant a s (/z 2 ). Recent 
investigations have shown that near threshold there can be large logarithms 
in the perturbation expansion which have to be resummed to make more 
reliable theoretical predictions. The application of these ideas to fixed-target 
Drell-Yan production has been the subject of many papers over the past few 
years ||. The same ideas on resummation were applied to the calculation of 
the top-quark cross section at the Fermilab Tevatron in and M. What 
is relevant in these reactions is the existence of a class of logarithms of the 
type (ln(l — z)) 1 / (1 — z), where i is the order of the perturbation expansion, 
and where one must integrate over the variable z up to a limit z — 1. These 
terms are not actually singular at z = 1 due to the presence of terms in 
5(1 — z). However the remainder can be quite large. In general one writes 
such terms as "plus" distributions, which are then convoluted with regular 
test functions (the parton densities). 

In this paper we examine the production of 6-quarks in a situation where 
the presence of these large logarithms is of importance, namely in a fixed- 
target experiment to be performed in the HERA ring at DESY. This actual 
experiment has the name HERA-B || ^| and involves colliding the circulat- 
ing proton beam against a stationary copper wire in the beam pipe. The 
nominal beam energy of the protons is 820 GeV, so that the square root of 
the center-of-mass (cm.) energy is \J~S = 39.2 GeV. Taking the 6-quark mass 
as rrib = 4.75 GeV/ c 2 then the ratio of mj/ y/S ~ 1/8. If we choose the renor- 
malization scale in the running coupling constant as m;, then a s (ml) ~ 0.2 
so a s (m 2 ) \xi(\fS / m^) m 0.4. This number is small enough that we expect a 
reasonably convergent perturbation series. 

In perturbation theory with a hard scale we can use the standard expres- 
sion for the order-by order cross section in QCD, namely 



a(S,m 2 ) = / dx x ^dx2^f i (x 1 ,fi 2 )fj(x2,f^)ffij(s=xix 2 S,m 2 ,fi 2 ), 



where the fi(x,fi 2 ) are the parton densities at the factorization scale /i 2 and 
the <7y are the partonic cross sections. The numerical results for the hadronic 
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cross sections depend on the choice of the parton densities, which involves 
the mass factorization scale fi 2 ; the choice of the running coupling constant, 
which involves the renormalization scale (also normally chosen to be // 2 ); and 
the choice for the actual mass of the 6-quark. In lowest order or Born ap- 
proximation the actual numbers for the cross section show a large sensitivity 
to these parameters. In section 3 we will show plots of the production cross 
section in leading order (LO), i.e. 0{a 2 s ) 1 and next-to-leading order (NLO), 
i.e. 0(al). The NLO results follow from the work of the two groups and 
[§, ||. However even including the NLO corrections does not completely fix 
the cross section. The sensitivity to our lack of knowledge of even higher 
terms in the QCD expansion is usually demonstrated by varying the scale 
choice up and down by factors of two. In general it is impossible to make 
more precise predictions given the absence of a calculation in next-to-next- 
to- leading order (NNLO). However in specific kinematical regions we can do 
so. 

The threshold region is one of these regions. In this region one finds 
that there are large logarithms which arise from an imperfect cancellation 
of the soft-plus-virtual (S+V) terms in the perturbation expansion. These 
logarithms are exactly of the same type mentioned above. We will see in 
section 3 that the gluon-gluon channel is the dominant channel for the pro- 
duction of 6-quarks near threshold in a fixed-target pp experiment. This is 
not the case for the production of the top quark at the Fermilab Tevatron, 
which is a proton-antiproton collider, and where the dominant channel is 
the quark-antiquark one. That was fortunate as the exponentiation of the 
soft-plus-virtual terms in || is on a much more solid footing in the qq chan- 
nel, due to all the past work which has been done on the Drell-Yan reaction 
[§]. Since the gluon-gluon channel is now the most important one we are 
forced to reexamine all "large" corrections near threshold, including both 
Coulomb-like and large constant terms. We will do that in section 2 where 
we will present all the relevant formulae at the partonic level and will also 
discuss the exponentiation of these terms. In addition we will present sub- 
leading S+V terms and discuss their contribution to the total S+V cross 
section. Section 3 contains the analysis of the hadron-hadron cross section 
which is relevant for the HERA-B experiment as well as for fixed-target pp 
experiments in general. We give results in LO, in NLO and after resumma- 
tion. Finally in Section 4 we give our conclusions and discuss where more 
work should be done in the future. 
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2 Results for parton-parton reactions 



The partonic processes that we examine are 

i{h)+j{k 2 )^Q(p 1 )+Q{p 2 ), (2.1) 

where i,j — g, g or i, j = q,q and Q, Q are heavy quarks (c, b, t). The square 
of the parton-parton cm. energy is s = (ki + A^) 2 - 

We begin with an analysis of heavy quark production in the qq channel. 
The Born cross section in this channel is given by 

a$/(s,m 2 ) = 2 -^al{^)K qq NC F \(3 (l + *f\ , (2.2) 

where Cf = {N 2 — 1)/(2N) is the Casimir invariant for the fundamental 
representation of SU(N), K qq - = N~ 2 is a color average factor, m is the heavy 



quark mass, /x denotes the renormalization scale, and (3 = y 1 — Am? / s. Also 
N = 3 for the 577(3) color group in QCD. The threshold behavior (s — > 4m 2 ) 
of this expression is given by 

m ") = ™ 2 s ^ 2 )K qq -NC F ±{3. (2.3) 

Complete analytic results are not available for the NLO cross section as some 
integrals are too complicated to do by hand. However in analytic results 
are given for the soft-plus-virtual contributions to the cross section, and 
for the approximation to the cross section near threshold. Simple formulae 
which yield reasonable approximations to the exact 0(a^) results have been 



constructed in |T0|. From these results one can derive that the Coulomb 



terms to first order in the qq channel are given by 

„<f>(W) = 4W)^ (C - (2.4) 



in the MS scheme, where Ca = N is the Casimir invariant for the ad- 
joint representation of SU(N). These terms are distinguished by their typ- 
ical behaviour near threshold which, after multiplication by the Born 
cross section, yield finite cross sections at threshold in NLO. We note that 
Cf — Ca/2 = —1/6 is a negative quantity for 577(3) and that the first-order 
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Coulomb correction is negative (the interaction is repulsive). From the orig- 
inal work of Schwinger [11] we know that the Coulomb terms exponentiate. 
Therefore we resum them by writing 



(7r 2 ), res / 
qq [ 



(0)/ 2\ 

<jqq J (s,m ) exp 



7r«s(/i s 
2(3 



Ca 
2 



(2.5) 



Since the exponent is negative, the exponentiation of the Coulomb terms 
actually supresses the cross section, and the resummed result goes to zero at 
threshold. In a previous treatment of threshold effects [E| the exponentiation 
was done in a different way, namely by writing 



where 



(7r 2 ),res/ 2\ (0)/ 2\ 

a g q {s,m ) = ay(s,m ) 
na s (fi 2 



(8) 



X, 



(«) 







c F - 



exp X(s) - 1 

C A \ _ 1 nets 
~ 6 (3 



(2.6) 



(2.7) 



The subscript in X(s) indicates that in the qq channel, where the process 
has to go via an intermediate s-channel gluon, the heavy quark pairs are 
produced exclusively in the color octet state. We have checked that the 
difference between the two methods of resummation is numerically negligible. 

In the DIS scheme in addition to the Coulomb terms we also have a large 
constant contribution so that the first order result near threshold is 



(t 2 )^ 2\ (0)/ 2\ 

a\q J (s,m ) = a q g'{s,m ) 



7ra s (/i 2 



c F - 



a s (jj? 



-C f 



7T 



9 

2 + y 



(2.8) 

We know that large constants exponentiate in the Drell-Yan reaction B so 
we assume that the same holds here and write the result in the DIS scheme 



as 



(tt 2 ), res 
T qq 



(s, m 2 ) 



(0)/ 2\ 

q q'{s,m )exp 



x exp 



7ra s (/i 2 
2(3 
'9 
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2 + T 



(2.9) 



This expression also goes to zero at threshold. We have included the constant 
terms to see their effect at larger values of (3. 
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Since the total parton-parton cross sections only depend on the variables 
s and m 2 they can be expressed in terms of scaling functions as follows 

oo 

<Tij(s, m 2 ) = ^^ fc) (s,m 2 ) 

k=0 

n; 2/,,2\ oo k 2 

= ^EC^WE/^fo)^ (2 ' 10) 

" L k=0 1=0 " L 

where we denote by the 0(a^ +2 ) contribution to the cross section. The 
scaling functions fij' l \rj) depend on the scaling variable rj = s/Am 2 — 1 = 
s(5 2 /Am 2 . 

In fig. 1 we plot fgq (t)) for k — 0, 1 for the exact and threshold ex- 
pressions (from ||) in the MS scheme. We also plot ^' Tes (r),a s ) which we 
define by 

(7r 2 ),res/ 2\ a s\^ ) /•(7r 2 ),res/ \ /o 1 i \ 

a\q h {S,m) = ^fqq ' 2.11 

We see that the threshold Born approximation is excellent for small rj and rea- 
sonable for the entire range of 77 shown. As expected the resummed Coulomb 
terms suppress the Born result throughout the entire range of i] and go to 
zero at threshold. We also note that the theshold first-order approximation 
is good only very near to threshold. In fig. 2 we plot the corresponding func- 
tions for the DIS scheme. Here the first-order corrections are smaller than in 
the MS scheme. Again the threshold first-order approximation is good only 
very close to threshold. The resummed result differs greatly from that in the 
MS scheme. Here the additional exponentiation of large constants cancels the 
negative contribution of the Coulomb terms and produces a large enhance- 
ment of the Born term in the region 0.1 < rj < 1 which, as we will show in the 
next section, is the most important region kinematically for the production 
of 6-quarks at HERA-B. For small rj, however, the total resummed result is 
dominated by the pure Coulomb terms and thus suppresses the Born term. 
Finally, we also show the resummed result where the only constant that we 
exponentiate is the 7r 2 /3 term in (2.9). In this case the enhancement of the 
Born term in our region of interest is much smaller. 

The analysis of the contributions to the gluon-gluon channel in NLO is 
much more complicated. First of all there are three Born diagrams each 
with a different color structure. Therefore only few terms near threshold 
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are proportional to the Born cross section. The exact Born term in the gg 
channel is 



Anai(fi 2 )K gg NC F - { C F 



1 

s 

Am 2 



\ S 

1 5 m 2 " 
3 + 3~ 



P + 





4m 2 \ 








s / 


m 4 \ , 


1 + /3" 


s 2 ) 111 


1-/3 


4m 4 

2 111 






1-/3 







(2.12) 



where X ss = (N 2 — 1) 2 is a color average factor. The threshold behavior 
(s — > 4m 2 ) of this expression is given by 



^t hres (^ 2 ) = ™ 2 s (v 2 )K gg -NC F [4C F - C A ](3. 



(2.13) 



Again, the complete NLO expression for the cross section in the gg channel 
is unavailable but analytic results are given for the S+V terms in || . These 
were used in [[10] to analyze the magnitude of the cross section near threshold. 
From the approximate expressions given in [10J one can extract the 7r 2 terms 
to first order in the gg channel. These are 



a£\s,m 2 ) = a^ 2 )NC K K gg ^- 



+ 32- 



m 



5 1 „ 9 rrr 

,4 



+a^CQEr>K gg - 



7T 



10- 
1 

"I 



m 



1, 1+/3 



16- 



m 



m 



m 



(2.14) 



where C K = (N 2 -1)/N = 2NC F C A - 4NC 2 F , and C QE d = (N 4 -1)/N 2 = 
—ACp + 4CaC f . These are not proportional to the Born term so that it is 
not clear how to resum them. The threshold behavior of (2.14) is given by 



( 2\ 

^93, thresh,™ 2 ) = a A s (n 2 )K, 



3/\.2\ 



99 4 s 



-NC 



K 



+ a 



QED 



(2.15) 
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which is proportional to the threshold Born term. Therefore the threshold 
approximation for the 7r 2 terms in NLO can be written as 



(0)+(tt 2 )/ 2 
°gg, thres I s ; 171 



a 



(0) 

gg, thres 



(s, m 2 ) 



1 + 



7ra s (/i 2 ) (~NC K /2 + C QED 



4(3 



(4CV - C A )NC F 



or, writing the color factors in terms of N, as 

7T« s (/i 2 ) N 2 + 2 



(0)+(tt 2 )/ 2\ _ (0) / 2 
^gg, thres { s , m ) ~ < J gg,thres\ S y m 



1 + 



(2.16) 



(2.17) 



4/3 N(N 2 - 2) 

The correction is positive indicating a Coulomb attraction in this channel. 
Then we proceed to resum these terms by the Schwinger method [11] 



a 



(tt 2 ), res 
gg, thres 



I 2\ (0) / 2\ 

(s,m ) = o-g a , thres (s,rn )- 



where 



X 



1 - exp(-X) 

7ra s (/i 2 ) iV 2 + 2 
2(3 N(N 2 - 2) 

42/3 v ; 



(2.18) 



(2.19) 



The reason that we have exponentiated in this way is that X is positive and 
tends to infinity when (3 — > 0. We expect that the Coulomb terms are only 
important very close to threshold. In |12| the Coulomb singlet (attractive 
interaction) and octet (repulsive interaction) contributions to the gg channel 
were exponentiated separately giving 



a 



is) 
gg 



-a 



(0). 



7 ^i-exp(-X (s) ) : 



X, 



4 note 



and 



(8) 



-a 



(0). 



(8) 



3 (3 
1 7ra s 



(8) 



6 (3 



(2.20) 



(2.21) 



7 99 exp(X {8) ) - 1 

respectively. Again we have checked that the difference between the two 
methods of exponentiation is not significant numerically. 

In fig. 3 we plot the scaling functions fgg(r]) with k — 0, 1 in the MS 
scheme for the exact and threshold expressions (from ||). We also show 
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fgg 2 ^ ,Tes which is defined in analogy to (2.11). We see that the Born and 
first-order threshold approximations are good only very close to threshold. 
The resummed result enhances the Born cross section and tends to a positive 
constant at threshold. The approximation is not good in the region 0.1 < 
77 < 1 so we turn now to a discussion of the important terms in this region. 

In m an approximation was given for the NLO soft-plus-virtual (S+V) 
contributions and the analogy with the Drell-Yan process was exploited to 
resum them to all orders of perturbation theory. The S+V approximation 
is adequate in the kinematical region of interest 0.1 < 77 < 1 for the qq 
channel, but not as good for the gg channel in the MS scheme. Therefore we 



reexamined the approximate formulae given in (Kj for the initial state gluon 
bremsstrahlung (ISGB) mechanism to see if there are subleading terms that 
will improve the S+V approximation. Let us see the structure of these terms. 
We are discussing partonic reactions of the type i{ki) + j(k 2 ) — > Q(pi) + 
Q{p 2 ) +g(k 3 ), and we introduce the kinematic variables t% = (k 2 —V2) 2 — m 2 , 
U\ = (fci — P2) 2 — in 2 , and s 4 = s + 1\ + u\. The variable S4 depends on the 
four-momentum of the extra partons emitted in the reaction. The first-order 
S+V result for the qq channel in the MS scheme is 



2 d 2 a^ (s,ti,ui) 
dt\du\ 



2C F 



71 



x 



r^l s 4 ,771" 

2 In -| + In — 

S4 \ m z jj, z 



+ 




0(s 4 -A) 



m 
J 2 



2C F s 4 



s 4 - A) 



In— b{s A ) 



2C 



where 



a 



qq 



(s,tx,ui) = vra s (/i )K qq NC F 



t\ + u\ 



+ 



2m 2 



(2.22) 



(2.23) 



s 

Here A is a small parameter used to allow us to distinguish between the soft 
(S4 < A) and the hard (s 4 > A) regions in phase space. The terms in the first 
pair of square brackets in (2.22) are the leading S+V terms given in |§ and 
those in the second pair of square brackets are subleading terms that we want 
to examine. In fig. 4 we plot the scaling functions /^'°' ) for the exact result, 
the leading S+V result, and the S+V result with both leading and subleading 
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terms. The leading S+V result is a reasonable approximation to the exact 
result in our region of interest 0.1 < rj < 1. The addition of Coulomb 
terms worsens the leading S+V result. We also see that when we include the 
subleading terms our approximation does not improve much in the region 
of interest. However, when we add both the first order Coulomb term and 
subleading terms to the leading S+V result we get a very good agreement 
with the exact result. Nevertheless, this is still not a major improvement 
over the simple leading S+V result. In the DIS scheme the analogous results 
are 



2 a °qq 
dt\du\ 



.3,*i,Ui) 



cr^(s,ti,Wi) a a {fi 




(2.24) 



In fig. 5 we plot the corresponding scaling functions. Here the addition 
of subleading terms worsens the leading S+V approximation. The addition 
of Coulomb terms and large constants enhances the first-order approximate 
results considerably. We also show, for comparison, the results of the addition 
of the Coulomb terms and the 7T 2 /3 constant term only to the approximate 
results. These last curves are the best fits to the exact result in the region 
0.1 < rj < 1. 

The resummation of the leading S+V terms has been given in ||. The 
result is 



, 2 dVSf(Mi,«i) 



dt\dui 



a 



qq 



s,h,ui) 



df(s4/m 2 , m 2 //i 2 



dsi 
A m 2 

+ /(— , — )*(*4) 
77T /J, z 



-6{s 4 - A) 



(2.25) 
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where 



/ s 4 


m 2 \ 




.C F _ 
A — a s 




In 2 4 






= exp 




\m 2 






7T 


\m z j 


m . 



2] V 



[sjm 2 ] 
V{l+r 1 ) 



exp(-?77 E J 



(2.26) 

Expressions for A, a s , rj, and 7# are given in [Q. As the NNLO cross section 
is not known exactly we do not how to resum the subleading terms. A 
reasonable guess would be (2.25) with the function / given now by 



./ 



s 4 ?7r 
m 2 ' ft 2 



Leading 



s 4 m 

m 2 ' /I 2 



exp 



Ca / 2\i S 4 

a 8 (fji ) In — 

7r m z 



(2.27) 



in the MS scheme, and 

,2\ 



/ 



s 4 m 
m 2 ' /i 2 



Leading 



s 4 m 

m 2 ' yU 2 



exp 



2tt 



3+41n2 - 



2C^ 



a s (n 2 ) In 



m 2 



(2.28) 

in the DIS scheme, where now we call /Leading the expression in (2.26). 

Now let us see the analogous results for the gg channel in the MS scheme. 
We have 



d 2 a^(s,t u 



dt\du\ 



2C, 



71 



1 

S4 



s 4 , m 
21n ^ + ln ^ 



+ 



3C, 



A , A 
^ + ln — 

m z m 

8C F 



0(s 4 - A) 

2\ " 



+5{s A ) I In 2 — + In — In — 

- r 1 i ■ i * ^ — 



-2C A + 8C F Vs 4 



I.0( S4 -A)+m— 5(s 4 ; 



(2.29) 



where 



2na 2 s {v 2 )K gg NC F 

t\ u\ 4m 2 s 
— + — A 



Cp — c , 



x 



m 2 s 



(2.30) 



Again, the terms in the first pair of square brackets in (2.29) are the leading 
S+V terms and those in the second pair of square brackets are subleading 
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terms. In fig. 6 we plot the scaling functions fgg for the exact result, the 
leading S+V result, and the S+V result with both leading and subleading 
terms. We note that the leading S+V approximate result is significantly 
smaller than the exact result and that the addition of subleading terms im- 
proves the approximation considerably. This is important since, as we will 
see in the next section, the gg channel is dominant for the production of 
6-quarks at HERA-B. Also the addition of Coulomb terms further improves 
the approximation. The resummation of the leading S+V terms for the gg 
channel has also been given in ||. The result is 



2 <Pal g s {s,t u ui) 



dt\du\ 



dfis^/m 2 , m 2 1 ' \i A 



+ f 



ds 4 
A 



(s 4 - A) 



m 



(2.31) 



where 



/ S 4 77T 

\m 2 ' jj? 





r c A 

2—a s 




In 2 il" 


I = exp 






7T 


\m z J 





r(i + r?) 



ex V (- VlE ). (2.32) 



Again, as the NNLO cross section is not known exactly we do not know how 
to resum the subleading terms. A reasonable guess would be (2.31) with the 
function / given now by 




/ S 4 77T \ 
/Leading! 2 , ^ I exp 



2Ca 3Ca — 8Cf , 2m s 4 
a s (/i )ln 



7T -2C A + 8C F 
where now we call /Leading the expression in (2.32). 



77T 



(2.33) 



3 Results for bottom quark production at 
fixed-target pp experiments and HERA-B 

In this section we discuss 6-quark production at HERA-B and also at fixed- 
target pp experiments in general, and we examine the effects of the various 
resummation procedures that were discussed in the previous section. Fol- 
lowing the notation in the total hadron-hadron cross section in order a k s 
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IS 

4 fc) (S,m 2 ) = E T dr^(r^ 2 )aS\rS,m 2 ^ 2 ), (3.1) 

where S is the square of the hadron-hadron cm. energy and i,j run over q, q 
and g. The parton flux $y(r, /i 2 ) is defined via 

^•(r,/i 2 )= ^.(a;,-^ 2 ), (3.2) 

Jt X X 

and ifjj is a product of the scale-dependent parton distribution functions 
fi(x, fi 2 ), where h stands for the hadron which is the source of the parton % 

H ij {x 1 ,x 2 ,v 2 )=f? 1 {x 1 ,v 2 )rf 2 {x 2 ,t J , 2 ). (3.3) 

The mass factorization scale /i is chosen to be identical with the renormal- 
ization scale in the running coupling constant. 

In the case of the all-order resummed expression the lower boundary in 
(3.1) has to be modified according to the condition sq < s — 2ms 1 / 2 , where 
so is defined below (see [|J). Resumming the soft gluon contributions to all 
orders we obtain 

a™ s (S,m 2 ) = £ f 1 dT$ i:i (T^ 2 )a tJ (TS,m 2 ^ 2 ), (3.4) 
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where is given in (3.24) of H] and 

[m+ (m 2 + s ) 1/2 ] 2 , Q _, 
r o = g , (3.5) 

with s = m 2 (nl/ jj 2 ) 3 ^ 2 (MS scheme) or s = m 2 (fj%/ii 2 ) (DIS scheme). Here 
fi is the non-perturbative parameter used in 0. It is used to cut off the 
resummation since the resummed corrections diverge for small fiQ. 

We now specialize to bottom quark production at HERA-B where v^S* = 
39.2 GeV. In the presentation of our results for the exact, approximate, and 
resummed hadronic cross sections we use the MRSD_' parametrization for 
the parton distributions Note that the hadronic results only involve 

partonic distribution functions at moderate and large x, where there is little 
difference between the various sets of parton densities. We have used the 



MRSD_' set 34 as given in PDFLIB [14[] in the DIS scheme with the number 
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of active light flavors n/ = 4 and the QCD scale A5 = 0.1559 GeV. We have 
used the two-loop corrected running coupling constant as given by PDFLIB. 

First, we discuss the NLO contributions to bottom quark production at 
HERA-B using the results in [7-9]. Except when explicitly stated otherwise 
we will take the factorization scale /1 = where m& is the 6-quark mass. 
Also, throughout the rest of this paper, we will use m and m& interchangeably. 
In fig. 7 we show the relative contributions of the qq channel in the DIS 
scheme and the gg channel in the MS scheme as a function of the bottom 
quark mass. We see that the gg contribution is the dominant one, lying 
between 70% and 80% of the total NLO cross section for the range of bottom 
mass values given. The qq contribution is smaller and makes up most of the 
remaining cross section. The relative contributions of the gq and the gq 
channels in the DIS scheme are negative and very small and they are also 
shown in the plot. The situation here is the reverse of what is known about 
top quark production at the Fermilab Tevatron where qq is the dominant 
channel with gg making up the remainder of the cross section, and gq and 
gq making an even smaller relative contribution than is the case for bottom 
quark production at HERA-B. The reason for this difference between top 
quark and bottom quark production is that the Tevatron is a pp collider 
while HERA-B is a fixed-target pp experiment. Thus, the parton densities 
involved are different and since sea quark densities are much smaller than 
valence quark densities, the qq contribution to the hadronic cross section 
diminishes for a fixed-target pp experiment relative to a pp collider for the 
same partonic cross section. 

In fig. 8 we show the K factors for the qq and gg channels and for 
their sum as a function of bottom quark mass. The K factor is defined by 
K = (cr(°) + c^ 1 ) |cxact)/o" ( ' ' ) , where is the Born term and | exact is the 
exact first order correction. We notice that the K factor is quite large for the 
gg channel, which means that higher order effects are more important for this 
channel than for qq. Since gg is the more important channel numerically, the 
K factor for the sum of the two channels is also quite large. We also show 
the K factor for the total which is slightly lower since we are also taking into 
account the negative contributions of the qg and qg channels. 

These large corrections come predominantly from the threshold region for 
bottom quark production where it has been shown that initial state gluon 
bremsstrahlung (ISGB) is responsible for the large corrections at NLO [[K] 



This can easily be seen in fig. 9 where the Born term and the 0(a s ) cross 
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section are plotted as a function of r] cut for the qq and gg channels, where 
i] = (s — 4m 2 ) /4m 2 is the variable into which we have incorporated the 
cut in our programs for the cross sections. As we increase 77 cut the cross 
sections increase. The cross sections rise sharply for values of i] cnt between 
0.1 and 1 and they reach a plateau at higher values of i] cnt indicating that 
the threshold region is very important and that the region where s >> 4m 2 
only makes a small contribution to the cross sections. This is the reason 
why we stressed in section 2 that our region of interest for comparison of the 
various approximations at the partonic level was 0.1 < rj < 1. In fig. 10 
we plot as a function of i] cut the Born term and the NLO cross section for 
the sum of the qq and gg channels and also the NLO cross section for the 
sum of all channels, including the small negative contribution of the qg and 
qg channels. We thus see that the qg and qg channels contribute a small 
negative contribution to the total exact NLO cross section. Note that in the 
last two figures as well as throughout the rest of this paper we are assuming 
that the bottom quark mass is m^ = 4.75 GeV/c 2 . 

Next, we discuss the scale dependence of our NLO results. In figs. 11 
and 12 we show the Born term, the exact first-order correction, and the total 
0(al) cross section as a function of the factorization scale for the qq and gg 
channels. We see that as the scale decreases, the Born cross section increases 
without bound but the exact first order correction decreases faster so that 
the NLO cross section peaks at a scale close to half the mass of the bottom 
quark and then decreases for smaller values of the scale. For the qq channel 
the NLO cross section is relatively flat. The situation is much worse for the 
gg channel, however, since the peak is very sharp and the scale dependence is 
much greater. Since the gg channel dominates, this large scale dependence is 
also reflected in the total cross section. Thus the variation in the NLO cross 
section for scales between m/2 and 2m is large. For comparison we note that 
the scale dependence for top quark production at the Fermilab Tevatron for 
m top — 175 GeV /c 2 is much smaller. 

In fig. 13 we plot the Born contribution for fx — m and the NLO cross 
section for fx = m/2, m, and 2m, as a function of the beam momentum for b- 
quark production at fixed-target pp experiments. The big width of the band 
reflects the large scale dependence that we discussed above. We see that 
the NLO cross section is almost twice as big as the Born term for the whole 
range of beam momenta that we are showing, and in particular for 820 GeV/ c 
which is the value of the beam momentum at HERA-B. The total NLO cross 
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section for 6-quark production at HERA-B is 28.8 nb for /x = ra/2; 9.6 nb 
for /x — m\ and 4.2 nb for /x = 2m. We also give the NLO results for the 
individual channels in fig. 14 for /x — m. 

In figs. 15 and 16 we examine the \x§ dependence of the resummed cross 
section for 6-quark production at HERA-B. We also show, for comparison, the 
/Xo dependence of + | app | app where we have imposed the same 

cut on the phase space of s 4 (s 4 > s ) as f° r the resummed cross section. 
Here | app and | app denote the approximate first and second order 
corrections, respectively, where only soft gluon contributions are taken into 
account. The effect of the resummation shows in the difference between the 
two curves. At small /xo, er rcs diverges, signalling the presence of the infrared 
renormalon. There is a region where the higher-order terms are numerically 
important, for instance 0.5 < /x < 1 in fig. 15. At high values of /x the 
two lines are practically the same. For the qq channel in the DIS scheme the 
resummation is successful in the sense that there is a relatively large region 
of /xo where resummation is well behaved before we encounter the divergence. 
For the gg channel, however, the situation is not as good. From these curves 
we choose what we think are reasonable values for /x . We choose /xo = 0.6 
GeV for the qq channel and /xo = 1.7 GeV for the gg channel. The value we 
chose for the gg channel is such that the resummed cross section is a little 
bit higher than the sum + | app | app . 

Using the values of /xo that we chose from the previous graphs, we proceed 
to plot the resummed cross section for 6-quark production at fixed-target pp 
experiments versus beam momentum. We present the results in fig. 17 for 
the qq and gg channels. We also show the results we have if in addition we 
resum subleading terms, and also if we resum both subleading terms and 
Coulomb and constant terms. For comparison the exact NLO results are 
shown as well. The resummed cross sections were calculated with the cut 
s 4 > s while no such cut was imposed on the NLO result. Since we know 
the exact 0{a z s ) result, we can make an even better estimate by calculating 
the perturbation theory improved cross sections defined by 

imp _ ^-res , _.(!) I _rr^ I C\ (\\ 

U H ~ U H ' U H lexact U H lapp i K -^) 

to exploit the fact that a$ | exact is known and | app is included in cr|f s . 
Then, in fig. 18 we plot the improved total cross section versus beam momen- 
tum (where we have also taken into account the small negative contributions 
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of the qg and qg channels) and, for comparison, the total exact NLO cross 
section for the three choices \x = m/2, m, and 2m. We also show the im- 
proved total cross section including resummation of subleading terms, and 
of both subleading terms and Coulomb and constant terms. The last curve 
lies above the NLO result for fi = m/2 for most values of the beam momenta 
shown (including the value at HERA-B) so that the effect of resummation 
exceeds the scale dependence of the NLO cross section. The improved total 
cross section for 6-quark production at HERA-B is 19.4 nb, if we resum lead- 
ing terms only; 23.8 nb, if we resum leading and subleading terms; and 31.4 
nb, if we resum leading, subleading, Coulomb, and constant terms. 

Finally, we present some results on the inclusive transverse momentum 
(p t ) and rapidity (Y) distributions of the bottom quark at HERA-B. Heavy 
quark differential distributions are known in NLO || Some of the rel- 
evant formulae for this part have been given already in []IE| , where the p t 



and Y distributions, including resummation, were presented for top-quark 
production at the Fermilab Tevatron. The heavy-quark inclusive differential 
distribution in pf is given by 

1^ 9l „ dr — J - — 2 , (3.7) 



dpt « 'J '4m»/5 dp 



with m\ = m 2 + p\. In the case of the all-order resummed expression the 
lower boundary in (3.7) has to be modified according to the condition s < 
s — 2m t s 1 / 2 . Resumming the soft gluon contributions to all orders we obtain 

daW(S,m 2 ,p 2 t ) fi zs dojjijS, m 2 ,p 2 ^ 2 ) 

"~ZW dT<bij[T,ii) — 2 , (3.8) 



dp 2 i? Jt * ' d P 



with daij/dp 2 given in (3.6) of JT^] and 



K + k + sq) 1 / 2 ] 2 

r o = g • (3.9) 

The corresponding formula to (3.7) for the heavy quark inclusive differ- 
ential distribution in Y is 

d<T%\S, m 2 ,y) /-i 2 da\f(TS,m 2 ,y^ 2 ) 

3T7 = 2^/ , 2 , dT$ij(T,p) — J - • 3.10 

dY ~j JAm 2 cosh 2 y/s dy 
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The all-order resummed differential distribution in Y is given by 

d* W (S,m>,Y) _ E ,1 dr$ij(T) ^ ) ^^ ) (3n) 



dY u ■' dy 



with daij/dy given in (3.9) of [|16| and 

[m cosh ?/ + (m 2 cosh 2 y + s ) 1//2 ] 2 
r = . (3.12) 

The hadronic heavy quark rapidity Y is related to the partonic heavy quark 
rapidity y by 

Y = y+\\n^. (3.13) 

We begin with the p t distributions. For these plots the mass factorization 
scale is not everywhere equal to m. We chose // = m in s , fu^s^/m 2 ,m 2 / ji 2 ) 
and a s , but /i = m t in the MRSD_' parton distribution functions and the 
running coupling constant a s (p). In fig. 19, we give the results for the qq 
channel in the DIS scheme. We plot the Born term da^ /dp t , the first order 
exact result da$ /dp t | exact, the first order approximation da% - /dpt | app , the 
second order approximation da H /dp t | app , and the resummed result da r ff /dp t 
for p, = 0.6 GeV. This is the same value for /x that was used above for the 
total cross section. We also show resummed results with the inclusion of sub- 
leading terms, and with both subleading terms and Coulomb and constant 
terms. If we decrease no the differential cross sections will increase. The re- 
summed distributions were calculated with the cut s 4 > s while no such cut 
was imposed on the phase space for the individual terms in the perturbation 
series. We continue with the results for the gg channel in the MS scheme. 
The corresponding plot is given in fig. 20. We note that the corrections in 
this channel are large. In fact the exact first-order correction is larger than 
the Born term and the approximate second-order correction is larger than 
the approximate first-order correction. In this case the value of /Jo has been 
chosen to be /xo = 1.7 GeV as above. We define the improved p t distribution 
by 

i i t | exact , I app • v'^ - -'- V 

dp t dp t dp t dp t 
In fig. 21 we plot the improved pt distributions for the sum of all channels, 
where we have included the small negative contributions of the qg and qg 
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channels. For comparison we also show the total exact NLO results for 
\i = m t /2, m t , and 2m t . The improved p t distributions are uniformly above 
the exact O(a^) results. We see that the effect of the resummation exceeds 
the uncertainty due to scale dependence. 

We finish with a discussion of the Y distributions. In this case we set 
the factorization mass scale equal to to everywhere. We begin with the 
qq channel. In fig. 22 we show the Born term da^/dY, the first order 
exact result da^ /dY | CX act, the first order approximation da^/dY | app , the 
second order approximation daff jdY | app , and the resummed result da T ^ s /dY 
for fi = 0.6 GeV. We also show resummed results with the inclusion of 
subleading terms, and with both subleading terms and Coulomb and constant 
terms. Again, the resummed distributions were calculated with the cut S4 > 
s while no such cut was imposed on the phase space for the individual terms 
in the perturbation series. We continue with the results for the gg channel in 
the MS scheme. The corresponding plot is given in fig. 23. Here, the value 
of Ho is fiQ = 1.7 GeV. The corrections in this channel are large as was the 
case for the p t distributions. We define the improved Y disrtibution by 

do$> _ daW , dv$ , doff 
dY dY + dY |cxact dY |app • 1 j 

In fig. 24 we plot the improved Y distributions for the sum of all channels, 
where we have included the small negative contributions of the qg and qg 
channels. For comparison we also show the total exact NLO results for 
H = to/2, to, and 2m. The improved Y distributions are uniformly above 
the 0(af) results. Again, we see that the effect of the resummation exceeds 
the uncertainty due to scale dependence. 



4 Conclusions 

We have presented NLO and resummed results for the cross section and dif- 
ferential distributions for bottom quark production at HERA-B. Results for 
the cross section as a function of beam momentum have also been given for 
fixed-target pp experiments in general. It has been shown that the gg channel 
is dominant and that the threshold region gives the main contribution to the 
NLO cross section. Approximations for the soft gluon contributions in that 
region have been compared with the exact results. The resummation of the 
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leading S+V logarithms produces an enhancement of the NLO results. The 
leading S+V approximation is not very good in the gg channel in the MS 
scheme in the kinematic region that is important for bottom quark produc- 
tion at HERA-B. The addition of subleading S+V terms and Coulomb terms 
improves the approximation considerably. The resummation of these addi- 
tional terms further enhances the cross section. We must stress, however, 
that our formula for the exponentiation of subleading terms is not based on 
any rigorous analysis and more work in this area will have to be done in the 
future. 
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Fig. 1. The scaling functions in the MS scheme. Plotted are f q q 
(exact, upper solid line at large 77; threshold approximation, upper dotted line 
at large if), fqg'°^ (exact, lower solid line at large rj; threshold approximation, 
lower dotted line at large if), and fqq ^' res (dashed line). 

Fig. 2. Same as fig. 1 but now for the DIS scheme. Also shown is fqq ' ) ' rcs 
where the only constant that we exponentiate is the n 2 /3 term (dash-dotted 
line). 

Fig. 3. Same as fig. 1 but now for f^'°^> in the MS scheme. 

Fig. 4. The scaling functions fgl' ^ in the MS scheme. Plotted are 
the exact result (solid line), the leading S+V approximation (dotted line), 
the leading S+V approximation plus Coulomb terms (short-dashed line), 
the S+V approximation with leading plus subleading terms (long dashed 
line), and the S+V approximation with leading plus subleading terms plus 
Coulomb terms (dash-dotted line). 

Fig. 5. Same af fig. 4 but now for the DIS scheme. Also shown are the 
leading S+V approximation plus Coulomb terms and the 7r 2 /3 constant term 
only (lower short-dashed line), and the S+V approximation with leading plus 
subleading terms plus Coulomb terms and the 7r 2 /3 constant term only (lower 
dash-dotted line). 

Fig. 6. Same as fig. 4 but now for /^'°- ) in the MS scheme. 

Fig. 7. Fractional contributions of the gg (MS scheme, short-dashed 
line), qq (DIS scheme, long-dashed line), qg (DIS scheme, lower dotted line), 
and qg (DIS scheme, upper dotted line) channels to the total O(a^) 6-quark 
production cross section at HERA-B as a function of 6-quark mass. 

Fig. 8. The K factors as a function of 6-quark mass for 6-quark production 
at HERA-B for the gg channel(MS scheme, short-dashed line), the qq channel 
(DIS scheme, long-dashed line), the sum of the gg and qq channels (dotted 
line), and the sum of all channels (solid line). 

Fig. 9. Cross sections for 6-quark production at HERA-B versus r] cut with 
mb — 4.75 GeV/c 2 for the qq channel in the DIS scheme and the gg channel 
in the MS scheme. Plotted are the Born term (gg, upper solid line at high 
i] cut ; qq, lower solid line at high ?7 cut ) and the 0(a^) cross section (gg, upper 
dashed line; qq, lower dashed line). 

Fig. 10. Cross sections for 6-quark production at HERA-B versus r] cut 
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with rrib = 4.75 GeV/c 2 . Plotted are the total Born term (solid line), the 
total 0(a^) cross section (dashed line), and the 0(al) cross section for the 
sum of the qq and gg channels (dotted line). 

Fig. 11. The scale dependence of the cross section for 6-quark production 
at HERA-B with = 4.75 GeV/c 2 for the qq channel in the DIS scheme. 
Plotted are the Born term (solid line), the exact first-order correction (dotted 
line), and their sum (dashed line). 

Fig. 12. Same as fig. 11 but now for the gg channel in the MS scheme. 

Fig. 13. The total Born (dotted line) and O(a^) (fj, — to solid line, \i = 
to/2 upper dashed line, and = 2to lower dashed line) 6-quark production 
cross sections at fixed-target pp experiments versus beam momentum for 
m b = 4.75 GeV/c 2 . 

Fig. 14. Contributions of individual channels to the total O(af) 6-quark 
production cross section at fixed-target pp experiments versus beam momen- 
tum for TO;, = 4.75 GeV/c 2 . Plotted are the contributions of the gg (MS 
scheme, short-dashed line) and qq (DIS scheme, long-dashed line) channels, 
and the absolute value of the contributions of the qg (DIS scheme, upper 
dotted line) and qg (DIS scheme, lower dotted line) channels. 

Fig. 15. The /i dependence of the resummed cross section for 6-quark 
production at HERA-B with to& = 4.75 GeV/c 2 for the qq channel in the DIS 
scheme. Plotted are a r q ^ s (solid line) and the sum + | app | app 
(dotted line). 

Fig. 16. Same as fig. 15 but now for the gg channel in the MS scheme. 

Fig. 17. Resummed and NLO cross sections versus beam momentum for 
6-quark production at fixed-target pp experiments for mj, = 4.75 GeV/c 2 . 
Plotted are the resummed cross sections for the qq channel in the DIS 
scheme for /x = 0.6 GeV (leading terms only, short-dashed line; with sub- 
leading terms, lower short-dash-dotted line; with both subleading terms and 
Coulomb and constant terms, upper short-dash-dotted line) and for the gg 
channel in the MS scheme for fiQ = 1.7 GeV (leading terms only, long-dashed 
line; with subleading terms, lower long-dash-dotted line; with both sublead- 
ing terms and Coulomb terms, upper long-dash-dotted line); and the 0(a^) 
cross sections for the gg channel in the MS scheme and the qq channel in the 
DIS scheme (upper and lower solid lines, respectively). 

Fig. 18. Improved and NLO cross sections versus beam momentum 
for 6-quark production at fixed-target pp experiments for mi, = 4.75 GeV/c 2 . 
Plotted are the total improved cross section (leading terms only, short-dashed 
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line; with subleading terms, lower long-dashed line; with both subleading 
terms and Coulomb and constant terms, upper long-dashed line) and the 
total 0{a\) result (/i — m solid line, /i — m/2 upper dotted line, \i = 2m 
lower dotted line). 

Fig. 19. The bottom quark p t distributions daffi/dpt at HERA-B for the 
qq channel in the DIS scheme for mj = 4.75 GeV/c 2 . Plotted are da^/dpt 
(solid line), da^/dpt | CX act (dotted line), da^/dpt \ app (short-dashed line), 
da H jdpt |app (long-dashed line), and da r ^ s /dp t for /io = 0.6 GeV (leading 
terms only, short-dash-dotted line; with subleading terms, lower long-dash- 
dotted line; with both subleading terms and Coulomb and constant terms, 
upper long-dash-dotted line). 

Fig. 20. Same as fig. 19 but now for the gg channel in the MS scheme 
and with /x = 1.7 GeV. 

Fig. 21. The bottom quark p t distributions dan/dpt at HERA-B for 
the sum of all channels for m b = 4.75 GeV/c 2 . Plotted are da^/dp t + 
da^/dpt | exact (/i = m t solid line, /x = m t /2 upper dotted line, \i = 2m t 
lower dotted line) and da^ p /dp t (leading terms only, short-dashed line; with 
subleading terms, lower long-dashed line; with both subleading terms and 
Coulomb and constant terms, upper long-dashed line). 

Fig. 22. The bottom quark Y distributions da^j jdY at HERA-B for the 
qq channel in the DIS scheme for mj = 4.75 GeV/c 2 . Plotted are da^ jdY 
(solid line), daff/dY | cxact (dotted line), daff/dY | app (short-dashed line), 
dcrjp jdY | app (long-dashed line), and da^/dY for /i = 0.6 GeV (leading 
terms only, short-dash-dotted line; with subleading terms, lower long-dash- 
dotted line; with both subleading terms and Coulomb and constant terms, 
upper long-dash-dotted line). 

Fig. 23. Same as fig. 22 but now for the gg channel in the MS scheme 
and with /i = 1.7 GeV. 

Fig. 24. The bottom quark Y distributions do H jdY at HERA-B for 
the sum of all channels for m b = 4.75 GeV/c 2 . Plotted are da^/dY + 
da^/dY | C xact (a* — 171 solid line, /j, — m/2 upper dotted line, /i = 2m 
lower dotted line) and da^ p / dY (leading terms only, short-dashed line; with 
subleading terms, lower long-dashed line; with both subleading terms and 
Coulomb and constant terms, upper long-dashed line). 
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